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Abstract.- The current density for a freely evolving state without negative mo- 
mentum components can temporarily be negative. The "operational arrival time 
distribution", defined by the absorption rate of an ideal detector, is calculated for a 
model detector and compared with recently proposed distributions. Counterintuitive 
features of the backflow regime are discussed. 
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The arrival time of a particle at a spatial point is one of the classical concepts 
whose quantum counterpart is problematic, even for the free particle case consid- 
ered in this paper. Despite this there are experiments, most notably "time of flight 
experiments", that seemingly circumvent the theoretical objections and difficulties 
exemplified by Allcock's work [l] . Several researchers have tried in recent years to fill 
this gap between theory and practice by reexamining the subject using a variety of 
approaches p|-^4|; a recent review provides a brief summary of the methods applied 
and a discussion of open questions [[14]]. 

If all of the particles in an ensemble of freely moving classical particles have pos- 
itive momenta p then the distribution n[t(X)] of arrival times t(X) at the position 
x = X is just the particle flux J[x = X,t]. (One spatial dimension is always assumed 
in this paper, the arrival point X is taken as the origin x = unless indicated oth- 
erwise, and t(0) is written simply as t.) Not surprisingly, in many investigations the 
proposed distribution H(t) of arrival times for quantum particles is closely related 
to the probability current J(0,t). In both the quantum and classical cases J(0,t) 
is equal to dN + (t)/dt, where N + (t) is the probability of finding the particle to the 
right of x = at time t (N + (t) = J °° \(x\ip(t))\ 2 dx in the quantum case). Hence, 
it is tempting to extend the classical result for the arrival time distribution to the 
quantum regime. But this is not entirely satisfactory because J(0,t) can be negative 
for a freely evolving quantum state even if its momentum distribution is zero for 
all negative momentum components, thus invalidating J(0, t) as a probability distri- 
bution of arrival times even in the free particle case. Bracken and Melloy showed 



that the time interval over which J(0, t) < can be arbitrarily long but finite [15 



They also derived a least-upper-bound, estimated to be 0.04, for the time integral 
of |J(0, t)| over an interval of negative J(0,t). But the backflow effect is negligible 
quantitatively at asymptotic distances from the source or interaction region [|J. This 
in part explains why the arrival time is not particularly worrysome for the practi- 
tioner of time of flight or other arrival time measurement techniques. However, the 
fundamental difficulty remains and is worth exploring. It is also expected that recent 
developments in atomic and optical physics will make the backflow regime amenable 
to experimental study. 

In this letter we will model an idealized particle detection setting and investigate 
the operational arrival time distribution focusing on the backflow regime. We have 
in mind a scintillation screen or any other device where the detection depends on 
the passage from the initial channel to one or more final channels (associated with 
changes in chemical arrangements or internal states of the particle or the apparatus). 
It is assumed that the experiment is repeated with single particles many times, with 
the same initial conditions, and that the final number of detection counts for a given 
dt is proportional to the amount of norm of the initial channel that disappears in 
that time, —dt(dN(t)/dt), where N(t) = dx \{x\ip(t))\ 2 , and (x\ip(t)) represents 
the amplitude of the initial channel. In general, in the absence of backflow at x — 0, 
the flux J(0, t) at the front edge of the detector (conventionally located between 
and L) and —dN(t)/dt are close to each other, but the latter is slightly delayed 
with respect to the former because of the time it takes to absorb (i.e., to pass from 
the incident to the final channels) the part of the wave inside the detector. More 
precisely, the time averages evaluated with J(0, t) and —dN(t)/dt differ by the mean 
dwell time in the detector m @- The arrival time distribution is in this context 
defined operationally, it depends on the apparatus, and it is given by the absorption 
rate —dN(t)/dt (suitably normalized to account for any incident particles that do not 
reach the detector). In real detectors additional delays and signal broadening have 
to be considered because of the amplification of the microscopic signal, for example 
in a photomultiplier, but we shall ignore that stage of the process, which is highly 
dependent on the particular detection method or apparatus, to concentrate on the 
generic first microscopic step. 

The effect of such a measuring device on the incident channel can be characterized 
by complex reflection and transmission amplitudes, R(p) and T(p), that depend on 
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the momentum. Within the spirit of "optical models" a complex potential may be 
constructed subject to the constraint of generating functions R(p) and T(p) with 
specified properties. This phenomenological approach retains the basic features of 
the apparatus-particle system avoiding the explicit treatment of all the degrees of 
freedom involved. A good detector will absorb completely over a broad range A p of 
incident momenta, 

R(p)=T(p) = 0, (pin A,). (1) 

With perfect amplification, all particles with incident energies in the working range 
will then be detected. Another desirable feature of a "time of arrival detector" is that 
its spatial width L be small. Clearly, these ideal conditions are difficult to meet in 
actual detectors, and known complex potentials do not satisfy them exactly either. 
It is however rewarding to study theoretical models as close as possible to the ideal 
limit to ascertain what to expect of "ideal measurements" . 

Consider an initial (t = 0) wavefunction with no negative momentum components 
and negligible overlap with the complex potential region [0, L \. In the absence of this 



potential the freely evolving wavefunction can be written as |JT6 

(x\4Hn(t)) = r dp{x\p)e- i ^ 2mh \p\<j )in {<S)) , (2) 

J 

whereas, if the potential is present the wavefunction is given by 

(x\m) = rMx\P + )e- ip2t/{2mh) (p\M0)) , (3) 



where the \p + ) are scattering eigenstates of H associated with the incident plane 
waves \p), and \(f>i n ) is the freely evolving incoming asymptote to which the state 
\if)) tends before the collision. The validity and meaning of this equation are not 
trivial. The elements of scattering theory with complex potentials required to obtain 
(HD are described in the Appendix B. For a "perfect absorber" the reflection ampli- 
tudes R(p) are zero for those momenta for which (p|0m(O)) is nonnegligible, so for 
x < the plane wave and the eigenstate of H coincide, (x\p) = (x\p + ). Therefore, 
(x\if>(t)) = (x\4>i n (t)) for x < 0, namely, the wave function to the left of the model 
detector is unaltered by the presence of the perfectly absorbing complex potential. 
The corresponding fluxes for and if) are also equal up to the potential edge, 

J^(x,t) = J<p(x,t), x<0. (4) 

This means that it is possible for a perfect absorber to emit probability! This coun- 
terintuitive phenomenon happens whenever probability backflow occurs at x = in 
the absence of the complex potential, i.e. J^x = 0, t) < for some finite time 
range 5 t . Since the perfectly absorbing detector reproduces this negative flux re- 
gion exactly it must "give back" part of the probability that had entered [0, L] 
before 8 t : dN~{t)/dt = —J^(x = 0,t) = —J^x — 0,t) > for t in 5 t , where 
N~(t) = J ^ \{x\if)(t))\ 2 dx. There is, however, no paradox because the leftward 
probability flow out of the detector volume is only temporary and does not imply a 
permanent reflection for a fraction of the particles since eventually all particles are 
absorbed. Thus a perfect absorber defined according to (|l|) does not require that 
J(0, t) always be nonnegative at its front edge. Moreover, the fact that dN(t)/dt < 
at all times for a complex potential with a negative imaginary part is not necessarily 
inconsistent with dN~(t)/dt being positive for a finite range St of t. These surprising, 
non-classical results will be demonstrated for an explicit choice of state and complex 
potential. The complex potential model is constructed by means of a series of com- 
plex square barriers with negative imaginary parts, see Appendix B. This enables us 
to know the scattering eigenfunctions exactly, and to calculate easily the absorption 
rate for a given wavepacket by quadrature (also quantities such as the total absorp- 
tion, or the dwell time). The method allows for efficient absorption at moderate to 
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large momenta (where Eq. ([!]) can be essentially satisfied with absorptions higher 
than 99.9%) but not for very low momenta. The wave packet is selected accordingly 
without too low momenta, and as a consequence the amount of backflow that we 
shall study is far from the upper bound found by Bracken and Melloy fl5 |, but this is 
not essential for demonstrating the effect and for illustrating the general theoretical 
prediction of Eq. @. 

We may first consider an incident state that leads to analytical expressions 
for (x\ip(t)) and J(x,t), and that satisfies the stated restrictions of the arrival time 
theories of references 
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(5) 



with 5 2 + a > and xq << — <5. Note the absence of negative momentum components 
and the p 2 dependence as p — > 0. The Fourier transform can be compactly expressed 
using w— functions [lgfl, 



(x\m) = c 
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where k = po/h. The derivative with respect to x, and therefore the flux, are also 
analytical since dw(z)/dz = —2zw[z) + 2i/ix 1 / 2 . The momentum distribution of this 
wave packet is too close to zero energy to be efficiently absorbed, so we shall use 
instead a new state obtained from (|) by a "boost" at time t = 0, 

(p|^(0)) = (p-6|V(0)). (10) 

The flux for the new state can be related to the flux and probability density of the 
original one, 

«V(0,t) = J4-bm/t,t) + (b/m)\(x = -bt/m\^(t)}\ 2 . (11) 

The last term may be understood as the flux contribution due to the displacement 
of an observer with velocity —b/m. Since b > 0, it is always positive, so if the 
boost is strong enough the backflow region disappears. We have chosen a value of b 
small enough that the effect remains and large enough that essentially full absorption 
can be achieved (99.97%). The complex potential has been improved in successive 
optimizations until the prediction of Eq. (f|) has been numerically confirmed by the 
converged flux results. Figures 1 and 2 show the flux at x = (with and without 
complex potential), and —dN(t)/dt. The two fluxes, with and without complex 
potential, are indistinguishable on the scale of the figures. This means in particular 
that, during the backflow regime detailed in Figure 2, the potential region [0, L] 
returns probability to the left half space [— oo, 0] at exactly the rate required to 
maintain the same negative flux as the freely evolving state. Fig. 1 shows global 
agreement in shape between the flux and the absorption rate —dN/ dt, the latter being 
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slightly delayed. It is of interest to compare the operational arrival time distribution 
—dN(t) I dt with other theoretical proposals for the time of arrival distribution. Also 
shown in Figures 1 and 2 are Kijowski's distribution 



n K (t) 



(rah) 1 / 

and that derived from Bohm's causal theory p~9|j20[1 



^e- ip2t/2mh (p\tP(0))dp 



(12) 



_ i T2[ [The transformation ( [T~T]) 
re non linearity introduced by 



U B (t) = \J(0,t)\/ \J(0,t')\dt', (13) 
Jo 

both evaluated for if)'. The former distribution was originally derived by imposing a 
series of conditions consistent with the classical distribution fl?fl , and has been later 
studied, rederived or generalized by several authors _ 
does not hold when J is replaced by Uk because of t 

the square root in ([12])]. It arises naturally as the square of the overlap between 
the initial state (restricted to positive momenta) and the eigenstates of the "time of 
arrival operator" 

m (J. 1 A 
2 \ p p J 

[The different quantization in || gives the same result.] As for any other quantity 
obtained from a formal quantization procedure, its physical meaning and content is 
not immediately obvious and requires detailed examination and external justifica- 
tion, since in principle other quantizations with the same classical limit may also 
be constructed [jSUl4j. In the numerical example Kijowski's distribution U.K(t) is in 
overall excellent agreement with the flux (except at the fine scale needed to resolve 
the backflow region) and, up to the delay, with —dN(t)/dt. With the latter it has in 
common positivity at all times. The two quantities avoid the negative values of the 
backflow region smoothly while |J(0, t)\ has downward cusps at the zeroes of J(0,t). 
In Bohm's theory particle trajectories do not intersect each other so that only a single 
trajectory contributes to the current density J(x,t) at each space-time point (x,t). 
For times t before S t , when J(0,t) > 0, particles in the ensemble crossing x = 
do so only from left to right; for times t during S t , when J(0, t) < 0, some of these 
particles recross x = 0, moving from right to left this time - there are no particles 
crossing x = from left to right during St] for times t after 8 t , when J(0, t) > again, 
particles cross x = only from left to right. By continuity, at those instants of time 
t when J(0, t) = no particle in the ensemble is crossing x = from either direction, 
leading to the above mentionned cusps. From a classical point of view, the coun- 
terintuitive aspect of this picture of backflow, when applied to an ensemble of freely 
evolving particles, is the necessity for some of these "free" particles to (twice) come 
to rest and then reverse their direction of motion. However within Bohm's theory 
the particles are not truly free but are guided by the wave function ip(x, t) [|2TJj which 
itself undergoes relatively rapid temporal changes in the backflow region considered 
here. Kijowski's distribution also contains an obvious counterintuitive feature and 
is subject to an interpretational puzzle. In different publications ||T7|.|5|. [T0| . [T2l , which 



all lead to the arrival time distribution fll2j) for a freely evolving state containing 
no negative momentum components, it is claimed that under those circumstances 
particles arrive at x = only from the left. In the backflow regime this leads to the 
baffling conclusion that the probability of finding the particle to the left of x = is 
increasing with time in a time interval during which particles reach x = only from 
the left. According to Bohm's theory, the particles that contribute to the backflow 
region shown in Fig. 2 arrive at x = three times, twice from the left and once 
from the right, so IIb(£) is not a "first arrival" time distribution; according to the 
approaches leading to Kijowski's distribution (O) they arrive only from the left but 



5 



it is not obvious how this is consistent with dN~ (t) jdt > and whether or not one 
should regard ( |12|) as a distribution of first arrival times. Finally, the operational 
distribution —dN/dt > has the advantage of a clear physical content, and direct 
relation to an experimental setting, but it should not be overinterpreted, in particular 
since the present analysis has shown that a perfect absorber can emit probability 



Appendix A: Construction of complex absorbing 
potential 

Let us find a potential with support [0, L] that maximizes the absorption in a given 
momentum interval [pi,P2]- The method described here is more efficient and numer- 
ically robust than previous ones ||22|| . The proposed functional form is a series of 
equal length N complex square barriers with complex energies {Vj}, j = 1,2, ...,N. 
The real and imaginary values of Vj are found by minimizing (with the restriction 
Im(Vj) < 0), the sum of the survival probabilities at s values of p, {p a }, 

f{V x , V N ;pi, ...,p s ) = £ S(V U V N -p a ) , (15) 

a=l 

where SOp) = 1 — \T(p)\ 2 — \R(p)\ 2 , and the s points are evenly spaced in the absorption 
interval [pi,P2]- S and its gradient with respect to {Vj} are obtained by multiplication 
of known 2x2 transfer matrices so that the optimizations are very fast. For the 
application in the main text we have taken L = 0.01, N = 4, s = 49, p\ = 260, and 
p 2 = 740. 



Appendix B: Elements of complex potential scat- 
tering theory in one dimension 

For the formulation of scattering theory of complex non hermitian potentials it is nec- 
essary to consider scattering eigenstates of H, \p + ), and their biorthogonal partners 
\p*~), which are eigenstates of H\ 

^ s W ^^-gt ^- < 17) 

Moller operators that connect the actual state ip that evolves, respectively, with H 
or W , with the freely evolving incoming asymptotic states (to which if> tends before 

the collision) can be defined as Q + = J dp \p + )(p\ and Cl + = J dp \p + )(p\ respectively. 
They obey 



it s 

it 



= i op , (i{ 



= i op -a, (19) 

where A = J2j is given in terms of the bound states of H (\F/), and of H\ 

i^j). Thus the integral form of an arbitrary wave packet evolving with H without 
bound state component can be written as 

(a#(t)) = / dp(x\p + )e- iE v t/h (p+\4){0)) . (20) 
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Using the generalized isommetry relation (Pf) one may substitute (p^^KO)) by (p|0m(O)) 
in (p0|) . Moreover, if ip(0) does not significantly overlap with the potential and there 
are not negative momenta, 0j n (O) can also by substituted by -^(O) as shown in p3| . 
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FIGURE CAPTIONS 



Figure 1: J(0,t) for free motion (solid line), J(0, t) with the absorber, (short 
dashed line), n^(t) (dotted-dashed), and —dN/dt (long dashed line) for the wave 
packet ip' in fllOD , see also (||), with the following parameters (all in atomic units): 
a = 1.4, p — 1, Xq — —0.22, 5 = 0.007, b = 300. The first three curves are 
indistinguishable in this scale. 

Figure 2: J(0, t) for free motion (solid line), J{0,t) with the absorber (short 
dashed line), |J(0, t)\ (dotted line), Hxit) (dotted-dashed line), —dN(t)/dt (dashed 
line), and -dN(t)/dt\ 

y+r D (dashed line with squares). Same parameters as in Figure 
1. The first two curves are hardly distinguishable, td = 1.0515 x I0~ 5 au. 
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